Diamagnetic response due to localization in chains of connected mesoscopic rings 
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A dynamic response to a magnetic field of a chain of con- 
nected mesoscopic rings is considered. We show that the low 
frequency behavior corresponds to localization of the elec- 
trons along the chain and to diamagnetic dynamic currents 
inside rings. The magnetization density due to these currents 
does not vanish even in the limit of the infinitely long chain of 
strongly connected rings showing that this is a macroscopic ef- 
fect. Being of a dynamic origin the currents can be destroyed 
by inelastic scattering but we argue that the corresponding 
decay time can at low enough temperature be large compared 
to all other time scales of the system. 
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The recent measurement of the magnetization of an 
array of disconnected Cu-rings Q constituted the first 
detection of a thermodynamic mesoscopic effect and re- 
vealed its dependence on the statistical ensemble. Con- 
siderable effort has been taken in order to build a theory 
of the observed effect in the diffusive regime where the 
circumference of a ring L is larger than the elastic mean 
free path I. Averaging over disorder at fixed Fermi energy 
gives an amplitude for the persistent current which is ex- 
ponentially small in L/l 0, and it was suggested that a 
larger amplitude can be due to the fixed number of par- 
ticles in the rings However, averaging over impurities 
in the canonical ensemble can analytically only be done 
by a perturbation theory in fluctuations of the chemical 
potential ||. 

Alternatively, it was proposed that the averaged dy- 
namic response to the magnetic field is less sensitive to 
the statistical ensemble Levels in a closed ring as 
functions of the external magnetic field do not intersect 
each other and, before the averaging, the dynamic re- 
sponse must coincide in the limit lo — > 0, where oj is fre- 
quency, with the corresponding thermodynamic deriva- 
tive. Therefore, in this case, one can hope that calculat- 
ing the dynamic response with a fixed chemical potential 
at T — OK, as was done in |Q, one obtains a reasonable 
description of the thermodynamic currents. Both the 
methods give the same order of magnitude of the per- 
sistent current in the diffusive regime [|| and the same 
result in the quasi-ballistic regime p|. Although, pos- 
sibly, not being exact, the dynamic approach can serve 
as a reasonable approximation scheme for calculation of 
thermodynamic quantities. 

However, let us emphasize that calculation of the dy- 



namic response in the isolated rings, due to the level re- 
pulsion, does not bring any new physics with respect to 
thermodynamics and can serve as a calculational tool 
only. (Recently, it was argued |j) that in a different ex- 
perimental setup the dynamic response can give results 
completely different from the thermodynamic ones). 

Now, while level intersections are not possible in an iso- 
lated ring, they can occur in a chain of connected rings in 
the limit of L x — > oo, where L x is the length of the chain. 
This is due to the fact that all wave functions in a disor- 
dered one dimensional chain are localized and parts of the 
chain located far from each other do not interfere. This 
leads to an exponential suppression of the level repulsion. 
The resulting level intersections in the system cause the 
dynamic response now to be completely different from 
the thermodynamic one. In fact, in the infinite system 
one cannot expect any finite thermodynamic magnetiza- 
tion density except the Landau diamagnetism which is 
rather small and will not be considered here. 

In this Letter we consider the dynamic response of a 
chain of connected rings in the limit of infinite length 
L x . We assume that every single ring is in the diffusive 
regime. It will be shown that in this disordered system 
the magnetization caused by an applied magnetic field 
is proportional to the length L x of the chain and, de- 
pending on the tunneling amplitude between the rings, 
can be measurably large. It decays in time due to inelas- 
tic processes, only. This effect is intimately related to 
the phenomenon of localization in quasi-one dimensional 
disordered systems. The system considered below offers 
the unique possibility to study localization in a magneti- 
zation measurement of an isolated system rather than in 
a transport measurement of the system coupled to leads. 

To understand better the nature of the dynamic re- 
sponse we want to study in the diffusive regime, let us 
consider first the corresponding clean limit. The chain 
of N r rings is in the ballistic regime L < I an ideal dia- 
magnet and its magnetization due to an applied flux (f> 
through the rings is given by: 

M dyn = NrVK^/Airc 2 , K d = -ne 2 /m, (1) 

where W is a small time dependent part of the flux 
through the ring. In Eq. (1) m is the electron mass, e 
-the electron charge, n -the electron density, c - the light 
velocity and V -the volume of a single ring. This mag- 
netization does not correspond to the free energy mini- 
mum and due to inelastic processes decays to the smaller 
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thermodynamic value which can exceed the Landau dia- 
magnetism, however Q. 

Let us now consider a model of non-interacting elec- 
trons moving in a random potential. The impurity aver- 
aged response to an external magnetic field H 



H = H + H u cos(u)t), 



(2) 



is calculated using the Kubo formula in the form given 
in Ref. [Q. The frequency dependent magnetization M u 
is then obtained from the response function K(lo) as 



M u = N r VRe [K{u))(j) u ] /4ttc 2 



(3) 



It will be shown below that in the absence of inelastic 
processes the dissipative part of the magnetization van- 
ishes in the zero frequency limit but the real part of the 
response can remain finite. In the limit of small frequency 
the diagrammatic technique leads to divergencies and one 
should use the supersymmetry method ||. 

The responses both along the chain of rings with cross- 
section S = V/L, as well as in the azimuthal direction 
of the rings can be reduced to a calculation of functional 
integrals over supermatrices Q. Here, we want to cal- 
culate the response in azimuthal direction of an ensem- 
ble of chains of rings whose length exceeds the localiza- 
tion length L x » L c >> L. The dynamic response for 
the problem of persistent currents in disconnected meso- 
scopic rings was calculated in Ref. B using the zero di- 
mensional version of the cr-model Q . In the general case 
of an arbitrary hopping amplitude between the rings one 
has to add in the free energy functional a term describ- 
ing the coupling between the rings. Such a term can be 
written in the form analogous to the Josephson coupling 
in the theory of superconductivity Jl(],[ll| ■ Changing the 
probability of tunneling from ring to ring we can describe 
thus the crossover from the case of isolated rings to the 
homogeneously weakly disordered chain of rings. Then 
the free energy F [Q] is written in the form 



F [Q] = - J2 J v STr IQiQA + Y. F ° \Q* 



(4) 
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where i and j enumerate the rings in the chain, Dq is the 
classical diffusion coefficient and we remind that Qf = 1. 
The notation STr stands for the Supertrace. The defi- 
nition of the matrices t% and A as well as details of the 
supersymmetry method can be found in previous works 
§,§. The function F [Qi] in Eqs.(|,|) stands for elec- 
tron motion inside the rings, while the first term in Eq.(|4|) 
describes the coupling between the rings due to electron 



hopping. The parameter J is related to the single elec- 
tron hopping amplitude t as J = ir 2 V 2 v 2 t 2 . Only nearest 
neighbors are coupled so that J„ = J for the neighbors, 
and Jij = otherwise. The limit J = corresponds to 
the chain of disconnected rings. In the limit J>1 the 
model on the lattice Eq. (Q) becomes the continuous one 
dimensional cr-model 0. Then, the coupling constant J 
is related to the classical diffusion constant Dq as 



J = irvVD /8L 2 



(6) 



Due to the one- dimensionality of the model Eq. (Q) one 
can use the transfer matrix technique. Corresponding 
partial differential equations have been written in Refs. 
^],[| for the continuous ID a- model. For the model on 
the lattice, Eq. (Q) , the corresponding recurrence equation 
is an integral one Jl^,|ll[. Repeating the main steps of 
these works we reduce calculation of integrals over Qi for 
all sites i to one integral over Q 



) Q = / (...)* 2 (Q)exp(-F Q (Q))dQ. 



(7) 

The function ^> (Q) in Eq. (Q) satisfies the equation 

* (Q) = J exp{2JSTrQQ')exp (-F (q')) * (q') dQ . 

(8) 

The solution of the Eq.(ph depends on the vector poten- 
tial A entering F , Eq.(@). In principle, it can be solved, 
at least numerically, for arbitrary magnetic fields using 
the parameterization proposed recently Q , but we will 
present results only in the limits of zero and high enough 
magnetic fields. These limits correspond to the orthog- 
onal and unitary ensembles. In both cases one can omit 
the first term in Eq.(||) because in the unitary case the su- 
permatrix Q commutes with the matrix T3 since the time 
reversal symmetry is broken Q . Due to the symmetry of 
the free energy the function 4" (Q) depends only on one 
compact real integration variable A and two non compact 
integration variables Ai t 2 parameterizing the supermatrix 
Q S. Therefore, in the functional integral over Q one 
can integrate first over all other variables of the repre- 
sentation of Q introduced in Ref. ||, reducing it thus to 
an integral over A and Ai^- 

Calculations for arbitrary frequency are most simple 
for the unitary ensemble, since the time reversal symme- 
try is broken and there is only one non compact integra- 
tion variable Ai. Corresponding calculations are not very 
different from those performed in Refs. 0] and we obtain 
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(uj) — — IUJ<Jq{1 



~ j J™ exp [(*™/A - 5) (Ai - A)] * 2 



(9) 
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where K umt (uj) is the response along the circumference 
of the rings in the unitary ensemble, A = {vV)~~ is the 
mean level spacing in one ring in the chain, <jq = e 2 Tn/m 
is the Drude conductivity and 6 — > 0. There remains a 
finite response after ensemble averaging, if at low fre- 
quencies uj the second term in the brackets in Eq.(^|) is 
proportional to l/u>. In the limit uj — > solving Eq.(||) 
becomes more simple because the main contribution to 
the integral over Ai comes from Ai ~ 1/uj. Then, the 
function '5 (Q) depends only on one variable z = 2ui\i 
and we obtain for the response K (to) 



{oj -+ 0) = -a A eff (J) /tt, 



(10) 



where the effective mean level spacing A e ff is a non- 
trivial function of the coupling J between the rings 



P(J)=A e// (J) /A 



exp (-z) ^> 2 j (z) dz (11) 



This function gives the probability that an electron does 
not leave a ring forever and is known numerically for 
arbitrary J [0- Here, we have to remark that the order 
in which the limits were taken, first taking the length of 
the chain L x to infinity and then performing the zero- 
frequency limit means that our results apply for chains 
of finite length only if the frequency uj is not smaller than 
the mean level spacing of the chain Ao = l/vSL x . 

At high frequencies the second term in Eq. (|^) is 
small and one obtains the classical response function 
K(ui) — — iujc o . In this limit the conduction is classi- 
cal in all directions. We see from Eqs. ( [| [h]) that the 
characteristic frequency of the crossover from the quan- 
tum to the classical regime is of the order of A e / / . 

In the limit J = corresponding to disconnected rings, 
^ = 1 and A e // (0) = A. The magnetization of the chain 
of disconnected rings is then, using Eq. (H): 
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f (J = 0) = NrVKdrA^/A^c 2 , 



(12) 



We remind that for the isolated rings the possibility of 
the averaging with the fixed chemical potential is a cer- 
tain assumption while for J ^ one should consider the 
whole chain and in the limit L x — > oo the difference be- 
tween the canonical and grand canonical ensembles van- 
ishes (however, taking into account charging effects one 
can have effectively isolated rings at finite J). 

In the opposite limit J > 1 the function ^ (z) is the 
solution of the differential equation 



zd 2 v/dz 2 - m.m = 0. 



(13) 



Solving Eq.(13) and calculating the integral Eq. (|1] 
we find that there is still a finite probability P(J) = 
1/(96 J) that an electron stays in one ring. This is a direct 
consequence of the localization of the electrons along the 
chain. In the unitary ensemble the localization length 
L c was calculated in Ref. || and can be related with the 



help of Eq.(||) to the coupling J as L c = 32 J L where we 
included an additional factor 2 since the effective cross- 
section of the chain of rings is 25. Note, that with an 
effective length of a ring L/2 in the chain the inverse 
participation ratio becomes T = 2/3L c which is the well 
known result in the unitary ensemble. In this limit we 
obtain for the effective level spacing A e ff (J) 

A eff (J>l) = (3vSL c )- 1 . (14) 

We see from Eq. (|l4| ) that in the limit of large J the 
response K umt (0) Eq. @ looks as if the chain of N R 
rings consisted effectively of L X /3L C (J) rings. In this 
limit we obtain with Eq. ([|) that the magnetization does 
not depend on the disorder: 

M umt {J ^ 1} = ^v/^Mt) 2 ) K^/Attc 2 , (15) 

where Mt — k 2 F S/A-K 2 is the number of transverse chan- 
nels in a single ring. 

Decreasing the coupling J makes the localization 
length shorter which leads to a larger response. 

The response K Eq. (10) describes the dynamic mag- 
netization per unit length along the chain and therefore 
the total magnetization is proportional to the length of 
the chain. It has been mentioned that the corresponding 
thermodynamic response is small as well as the differ- 
ence between the canonical and grand canonical ensem- 
bles. To make a more explicit estimate, one can use the 
expansion in terms of fluctuations of the chemical poten- 
tial ||,[l4| and calculate the correlation functions enter- 
ing the expansion with the supersymmetry technique. It 
turns out that the thermodynamic magnetization density 
of the chains calculated in this approximation vanishes 
for macroscopically long chains |L3| . 

In order to estimate the decay time of the dynamic 
magnetization one has to specify inelastic processes. As 
we found out above, the dynamic magnetization is due 
to the localization along the chain. As soon as the lo- 
calization is lifted and electrons can travel through the 
chain, the magnetization should vanish. At temperatures 
T > A e ff and at frequencies uj smaller than t^ 1 (J) , 
where t (J) is a time due to inelastic processes, the dy- 
namic magnetization decays due to a d.c. conductiv- 
ity proportional to t^ 1 (J) whereas at higher frequencies 
t(J) has no influence on this magnetization. The sim- 
plest way to include the time t ( J) into formulae obtained 
above is to substitute uj — > uj + it^ 1 (J) in all expressions 
for K/uj, where K can stand for the responses in both 
the longitudinal and azimuthal directions. (K/lu is pro- 
portional to the product of two Green functions and it 
is this quantity that is calculated in the linear response 
theory). With this substitution the magnetization would 
decay with time t (J). Of course, this is a hypothesis and 
we suggest it to make a rough estimate only. As concerns 
the longitudinal conductivity in the regime of the local- 
ization it can be different from zero due to hopping be- 
tween the localization centers [n5[. Explicit calculations 
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for disordered chains |16|] at not too low temperatures 
corresponding to T > A e ff lead to the result t (J) = r p h 
where r p h is the electron-phonon scattering time. To the 
best of our knowledge nobody has performed analogous 
calculations for thick wires which are equivalent to the 
chain of the rings but we hope that this estimate can give 
a reasonable estimate also for the case considered here. 
Then, the inverse decay time is proportional to T 3 . Thus, 
at temperatures much below the Debye temperature ujd 
but above the effective level spacing, uod » T > A e ff, 
one can have l/t(J) < A e ff and it should be possible to 
observe the large diamagnetic response. 

At lower temperatures T <C A e / / the exponential Mott 
law, t(J) ~ (l/w D )exp(4(3A e// /T) 1 / 2 ), might serve as 
a lower limit for the decay time. Then, it increases ex- 
ponentially with decreasing the localization length. 

We note that charging effects might increase both the 
amplitude and the decay time of the dynamic magneti- 
zation. A detailed analysis of these effects will be the 
subject of further research. 

Analogous calculations can be performed also for the 
orthogonal ensemble corresponding to the zero static 
component H of the magnetic field Eq. (||) . In this case 
we obtain as in Ref. Q 

M°Jl h (0) = 0. (16) 

Changing the static component H of the magnetic field 
we can have a crossover between Eqs.(|l6|) and (|i"o|). The 
characteristic flux <f> c of this crossover is of the order of 
0o(A e// /£: c ) 1/2 , where (j) = 2nc/e and E c = tt 2 D /L 2 
is the Thouless energy, cj) is the flux quantum. The de- 
pendence of the response on the flux <f> is periodic with 
the period c/>o/2. 

Adding magnetic or spin-orbit impurities changes Eqs. 
( |l0| , |l6| ). The system with the magnetic impurities 
corresponds to the unitary ensemble. One can use as 
before Eq. ( |io| ) provided A e // (J) is substituted by 
A e ff (2 J) /2. If the magnetic impurities are absent spin- 
orbit ones lead to a different function A e ff (J). However, 
this difference is only numerical and does not change the 
sign of the response which is in all cases diamagnetic. 

In conclusion, we showed that a magnetic field applied 
perpendicular to a chain of disordered connected rings in- 
duces a macroscopic diamagnetic current. In the absence 
of inelastic scattering which can correspond to low tem- 
peratures this current can live for a very long time. The 
corresponding magnetization density remains finite even 
in the limit of a macroscopically long chain of perfectly 
connected rings. At the same time the longitudinal re- 
sponse along the chains shows a dielectric behavior usual 
for disordered wires. In fact, the shorter the localiza- 
tion length along the chain is, the larger is the current 
in the rings and the larger is the relaxation time to its 
vanishing thermodynamic value. We suggest to check 
our theoretical results by measuring the magnetization 



as a function of flux of chains of mesoscopic rings in the 
diffusive regime or, alternatively, of linear antidot lat- 
tices where the mean free path due to elastic scattering is 
smaller than the distance between neighboring antidots. 
This can yield a new contactless method of studying the 
localization. 
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